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Abstract 


We calculate the masses of the nucleon resonances [N+(1440), N-(1535), N- (1650), N+(1710), N+(1880), N-(1895), N+(2100)] 
using an adaptation of the method of finite energy QCD sum rules (FESR) introduced recently in the determination of the 
masses of the vector meson resonances. This variation of finite energy sum rules involves as usual integration over a closed 
contour consisting of a circle of a large radius R and a cut in the complex energy (squared) plane. A specific integration kernel 
is used, however, which minimizes the integral over the cut. We obtain remarkably stable results over a wide range R, where 
R is the radius of the integration contour. The sum rule predictions agree well with the experimental values. We make use of 


the nucleonic correlator discussed in our previous FESR calculation of the nucleon mass. 


Keywords: Sum Rules; QCD; Nucleon Resonance Masses 


Abbreviations: OPE: Operator Product Expansion; FESR: 
Finite Energy Sum Rules. 


Introduction 


QCD Sum Rules and Nucleon Resonance Masses 


QCD sum rules based on the Borel transformation of 
the operator product expansion (OPE) were applied to the 
calculation of low energy properties of hardons starting with 
the pioneering paper of Shifman, et al. [1]. The method was 
soon afterwards applied to the nucleon by Furnstahl RJ, et al. 
[2,3]. An alternative more general approach is based on Finite 
energy sum rules (FESR) [4,5]. Most recently was applied we 
applied our approach to the calculation of all vector meson 
masses [6]. We loffe BL, et al. [2] here extend our work to the 
calculation of the masses of the nucleon resonances. 


Nucleon Resonance Masses from QCD Sum Rules 


Consider the nucleon correlator 
M(t = 4°) =i] dre (0|Tm(x)n(0)]0) 
where 
n=é (u,Cyu, )s¥ 24 


The nucleon current proposed in Ioffe BL, et al [2]. In 
a previous work Nasrallah NF, et al. [7,8] we studied this 
correlator to calculate the mass of the nucleon. We used 
FESR with polynomial Kernels designed to eliminate the 
contribution of the nucleon continuum. This continuum 
consists mainly of the 


N* (1440), N* (1535),.N* (1650), N* (1710), N* (1880), (1895), N* (2100) 
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Here we shall use a slightly modified integration kernel to 
calculate the masses of nucleon recurrences beginning with 
the Roper resonance N +(1440) (Equations 1-8). 


The amplitude IT(t) can be decomposed as 


T(t) = gl, (t)+T1,¢ (1) 


We shall work with IT, (¢). At low energies 


we (2) 


t—my 


II, (t) =—-A’m, 


where { is the coupling of the current to the nucleon. 
At high energies 


(27) 0 aca (t) ~ Bt In(-1) +224 (3) 
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The correlator and the parameters were given in Nasrallah 
NE et al. [7,8] where we calculated the nucleon mass proper 
and presented a detailed analysis 


B=ag (aa) 3.4) =-0.908GeV (4) 


2 
B, = ~~ (qa)(a,6") = 0.034Gev" (5) 


7 


B, =~(2z)° = (aay a, =0.083GeV® (6) 


a, 


where we use for a, = =0.10, the quark condensate 


a 


(aq) =-0.02GeV; and the gluon condensate 


a,G? =0.013GeV’" 


The basis of FESR is Cauchy’s theorem applied to the contour 
C of Figure 1 which implies 
lf I 
—A’m,P(m;,)=— | dtP(t)ImI1™ (t)+—@ dt P(t)T1(t) (7 
P(m;,) = (t) (t) ni |J=R (t) (t) (7) 


cut 


On the circle of large radius R the exact correlator 
T(t) can be replaced by its QCD expression 11° (r) . 


This is the usual definition of duality. The QCD correlator 
involves a perturbative part and a non-perturbative part 
which we parametrise by condensates [8]. We do not 
consider genuine duality violations. 


Nasrallah NF and Schilcher K. Nucleon Resonance Masses from QCD Sum Rules. Open J of Astro 


2024, 2(1): 000108. 


Figure 1: Integration Contour of FESR. 
see. Se 


Re(t) 


Where P(t) is an entire function, e. g. a polynomial. Over 
the circle of large radius R the correlator II (t) has been 
replaced by its QCD expression. The principal unknown in 
Equation 8 is the integral over the cut, i.e. over the higher 
nucleon resonances with mass my, <R. To minimize this 
integral (before neglecting it), a judicious choice of the 
weight-function P(t) has to be made. With the classic choice 
P(t)= exp(-1/M;) the Borel variable M@,; cannot be chosen 
too large because it would minimize the contribution of the 
nucleon. Also M; cannot be too small because the unknown 
condensates in Equation 3 would explode. It was hoped in 1 
that a region of stability at an intermediate /, 9 can be found. 
This can be shown to be not the case 3. In our FESR approach 


Copyright© Nasrallah NF and Schilcher K. 


Nasrallah NF, et al. [7] we take P (t) to be a polynomial 


= Sor" (8) 


n=0 


Itis clear that the order n, ,. cannot be chosen arbitrarily 
high because of the contribution of unknown _ higher 
condensates. To test the method we start with the nucleon 
proper (Equations 9-19). We choose 


P(t,R) ' wa : *| l-at-a,t? (9) 


where a,R= + 


Equation 7 then becomes 
(27) Amy P(my, ) = 


(27) A? my P (my ) = 


—B,I,(R)- By +A, (0) 
—B,1,(R)-B,+a4,B,+A, 11) 


where 7, = [att P(t,R), i; = [dtt? P(t,R) (12) 


and 


apr apr 
=-(2n)'[ deP(s)ImI1, (1), A,=-(27)'f dee P(0)imm, (1) 13) 
P (t) is designed to minimize the contribution of the 
continuum in Equation 7, i. e. to minimize the contribution of 
A, + A, so as to neglect it then. The nucleon mass is obtained 
from the ratio (Equations 10/ 11). 


There is a wide region of stability R for /,(R) and /,(R) at 
around R = 3 GeV’ Using the sum rule Equation 7 


(27) A?m,,P (m;,)=—B,I, (R)- B,—A, (14) 


(27)' A?my,P (m;,)=—B,1,(R)- B, +4,B,—A, (15) 
where 
A,(t)=-(2)' f* deeP(t)imMl, (1) and A, (t)=-(27)' [* dee P(e), (t) (16) 
The delta’s denote the contribution of the nucleon 
continuum, they receive no contribution from the nucleon 
pole [9]. The aim is to choose P(t) which minimizes A, and 
A, so as to neglect them and obtain m, from Equations (1/2). 
This was done in 7 from the choice 


P(t)=1-at—a,t* =1-0.807t+0.16t? (17) 


3GeV? 


This choice minimizes the integral Ves _dt |P(t ) 


For R = 2.5 GeV’ the result was 
my = 0.830+0.05GeV (18) 
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Note that to obtain the physical mass one has to add the 
contribution of the X-term. Using higher moments we get we 
get information on the higher resonances. 


We Nasrallah NF, et al. [7,8] next proceed with a different 
choice for P (t) which eliminates A, , as well. 


P\(,R)= pL -2) (19) 


which vanishes at the mass of the m, where m, =m, 
(1440) is the mass of the First exited state (the Roper) where 
we expect the bulk of the contribution to the continuum to 
come from. R will be determined by stability considerations 
(Equations 20-28). 


The expression for m, is 


1 (R) + 
Mm, = 8; : (20) 
1(R)+3h a, 


where 


I,=['dt? P(t,R), T=] dteP (t,R) (21) 


: ; ; 1 1 ‘ 1 
P(t)=1-at-at’, a = +—|la,= (22) 
(t) 1 2 1 [ +) 2 mR 


1 


With my as an input, Equation 20 determine m, asa function 
of R: 


R=2.1GeV? m,=1.57GeV 
R=24GeV? m,=146GeV 
R=3.0GeV? m,=149GeV 


We Nasrallah NF et al. [7,8] combine and quote 


m, =1.46+0.03GeV 


The procedure can be repeated using the general formula 


cde (MPP, (tR) 
ee att Pltk) 


eh aha t t 
wit (AR): nf +l *| (23) 


R being the maximum at which the maximum of 
ee Ht hae (t, R) takes place (which is also close to the 
maximum of ate’ Pleh) 
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We Nasrallah NF, et al. [7,8] list our results below 


m, = my (1535)=147GeV, R=3.2GeV" 

m, = my (1650) =1.57GeV, R=3.5GeV* 

m, = my (1710) =1.74GeV, R=4.0GeV" 
mM, , = My (1880,1895) =1.80GeV, R=5.3GeV" 


m, = my (2100)=1.98GeV, R=6.4GeV" 


The errors of the predicted masses of resonances are about 
+10%. 


Heavy Baryons 


The A,,A,,4,,=, and their parity doublets can be treated by 
the same sum rule method using an integration kernel of the 


form [10] 
(um?) =(1- 1-2) (24) 


where m’ is mass of the first exited state or parity 
doublet and R is the radius of the circle of Figure (1). 


We Nasrallah NF, et al. [7,8] first determine the mass of the 


ground state baryon fora knowledge of m’ and then determine 
m’ independently and self-consistently. We consider 


A, (2ss0s)[a.4 esss(od | 


. (sco) 05 J (orn) a.t | 


In the correlator of Equation 7 we use the simplest current 
for Ay 


We define as above 
; 1 
JO = V2 [abe (ui Cyd, —d7Cyzu,)Q.) (25) 
I1(t) = gA(t)+ Br (26) 


On the hadronic side we make explicit the contribution of 
the positive and negative parity baryons, 
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where 4. the coupling of the current Equation 25 to the 
corresponding states (Equations 29-36). 


The QCD spectral functions have been obtained in Nasrallah 
NF et al. [7,8], Ayala C, et al. [9] and Krasnikov V, et al. [10] 
and are conveniently expressed in Gorishnii SG, et al. [11], 
Kataev AL [12], Dominguez CA [13], Zyla PA, et al. [14], Wang 
ZG [15] and Zhao ZX, et al. [16]. 


We define 3m 
A — 
pa(i)= 22-0) 
with ; ; 
o()=[garstt-ay |S -t) «OS 22) eo 
mo My Xx 6 mM, t 


4 
6 mM 


1 
Consider now the integral Flot APH) where C is the 
contour of Figure 1 in the complex t-plane. 
We Nasrallah NF, et al. [7,8], Wang ZG [15] claim that our 
choice of the damping kernel P (t, R) essentially eliminates 


all hadronic contributions except that of the ground state so 
that 


22P(m?)= |" dtP(t,R)p,(t)=1, (30) 

mo 

Another equation is obtained from the First moment 
A2m?P(m?) =f" dttP(t,R)p,(t)=1, 


R is chosen in the stability region of the integrals J, and J, 
for A,. These integrals turn out to be very at functions of t 


which obtain a maximum between 4.5? and 5m: . The result is 
my, =2.10GeV (exp.2.23GeV) (31) 
my, =5.92GeV (exp.5.62GeV) (32) 

Using B(t)P (t, R) yields practically identical results. 

We now proceed to calculate m, independently. Because b 


and c quarks are heavy T1(0) is given by its QCD expression 
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a 
cy 


which a an additional sum rule 
i R at 
-P(ni)=[" P(.R)ay()=1, C9) 
m, mg. 


This implies 


from which we can determine m.. The result is 
m,. =(2.5140.03 )GeV (exp.2.58GeV ) (34) 
m, =(5.87+0.14)GeV (exp.5.91GeV) (35) 


The errors are estimated by varying R by 10%. 
Finally for the & baryon =, ~~ a Ey ae ue 
5797 \ 2 2 “(5935 ){ 2 2 


the function a(t) in Equation 29 picks up an additional term 
proportional to the strange quark mass m, 


3 ms 


2 Ts 2H 4 
Aa(t)="" ait ah : | 
Q 


The calculation proceeds as before, giving 


mz =2.05GeV (exp.2.468GeV ) 
mz, =5.40GeV (exp.5.79GeV ) 
m, =2.46GeV (exp.2.58GeV ) 


m. =5.87GeV (exp.5.935GeV) 


Conclusion 


We have calculated the masses of the baryon recurrences 
with a new variant of QCD finite energy sum rules. The only 
free parameter of the sum rules, the radius of the circle in 
the complex t-plane, is fixed by the requirement of stability. 
The method works well for all similar systems such as the 
vector resonances. The main source of error is the zero 
width approximation for the resonances. We have estimated 
this error by allowing the radius entering the sum rule to 
vary by +10%. Order «a, corrections are included, order 

s are calculated and found to be negligible. The sum rule 
predictions are compared with the experimental numbers 
and agreement within the expected accuracy is found. It can 
be concluded that QCD is applicable to single resonances and 
their recurrences. 
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